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Abstract. The Dirac Hamiltonian with relativistic scalar and vector harmonic oscillator po- 
tentials has been solved analytically in two limits. One is the spin limit for which spin is 
an invariant symmetry of the the Dirac Hamiltonian and the other is the pseudo-spin limit for 
which pseudo-spin is an invariant symmetry of the the Dirac Hamiltonian. The spin limit occurs 
when the scalar potential is equal to the vector potential plus a constant, and the pseudospin 
limit occurs when the scalar potential is equal in magnitude but opposite in sign to the vector 
potential plus a constant. Like the non-relativistic harmonic oscillator, each of these limits has 
r^ ■ a higher symmetry. For example, for the spherically symmetric oscillator, these limits have a 

r^ U(3) and pseudo-U(3) symmetry respectively. We shall discuss the eigenfunctions and eigen- 

values of these two limits and derive the relativistic generators for the U(3) and pseudo-U(3) 
symmetry. We also argue, that, if an anti-nucleon can be bound in a nucleus, the spectrum will 
>._^ ^ have approximate spin and U(3) symmetry. 
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(<— ^ . 1. Introduction 

C^ I Pseudopsin symmetry is approximately conserved in nuclei [1]. Pseudospin doublets have single 

nucleon quantum numbers [n,£,i + ■^] and [n — l,i + 2,i + ^] where n,£ are the radial and 
orbital quantum numbers, respectively, and the last quantum number in the brackets is the 
total angular momentum, j. For example, in the usual notation, (IsijOds), (IpsjO/s), etc, are 
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Vh ' pseudospin doublets. These doublets are approximately degenerate in energy. Some beautiful 



examples have been shown in this conference; for example, in the talk by O. Sorlin [2]. Most 
of the properties of nuclei are explained by the non-relativistic shell model [3]. However, to 
explain pseudospin symmetry we need the relativistic Dirac Hamiltonian which has two types 
of mean fields, a Lorentz scalar Vs{r) and Lorentz vector Vv{r). Pseudospin symmetry occurs 
when the sum of the vector and scalar potentials in the Dirac Hamiltonian is a constant , Vs{r) 
+ ^v(^) = Cps [1]. Hence pseudospin symmetry is a relativistic symmetry. The approximate 
equality in magnitude of the vector and scalar fields in nuclei and their opposite sign have been 
confirmed in relativisitc mean field theories [T] and in QCD sum rules [5]. 

On the other hand hadrons [IJ and anti-nucleons in a nuclear environment have spin symmetry 
[6]. These are relativistic systems and normally, in such systems, we would expect large spin- 
orbit splittings. However spin symmetry occurs when the difference of the vector and scalar 
potentials in the Dirac Hamiltonian is a constant ^, Vs{r) - Vvir) = Cg, even though these 
systems are very relativistic. 



The non relativistic harmonic oscillator has been useful in the shell model of nuclei. However, 
the relativistic mean field has been important for pseudospin in nuclei and for mesons with one 
heavy quark [8] and for an anti-nucleon in a nuclear environment [9]. For these reasons the 
relativistic harmonic oscillator has been solved in the tri-axial, axially deformed and spherical 
limit [iQj. Like the non-relativistic harmonic oscillator, each of these limits has a higher 
symmetry. For example, for the spherically symmetric oscillator, these limits have a U(3) and 
pseudo-U(3) symmetry respectively. We shall discuss the eigenfunctions and eigenvalues of these 
two limits and derive the relativistic generators for the U(3) and pseudo-U(3) symmetry. 

2. The Spherically Symmetric Dirac Hamiltonian with Spin Symmetry 

The Dirac Hamiltonian for a spherical harmonic oscillator with spin symmetry is 



H = a-p + (3M + {l + /3)V{r), (1) 

where a, (3 are the Dirac matrices, p is the momentum, M is the mass, r is the radial coordinate, 
r its magnitude, and the velocity of light is set equal to unity, c =1. The generators for the spin 
SU(2) algebra and the orbital angular momentum SU(2) algebra, S, L, which commute with the 
Dirac Hamiltonian with any potential V{r), [H,S] = [H,L] = 0, are given by [11] 



S 



s 

UpsUp 



L 



i 

UpiUp 



(2) 



where s = a/2 are the usual spin generators, a the Pauli matrices, £ = iZ^^H^ 
is the helicity unitary operator introduced in |12j . 
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2.1. The Dirac Hamiltonian with Spin Symmetry and Harmonic Oscillator Potential 



With the harmonic oscillator potential V{r) 



Moj^ 



r the eigenvalue equation is ^0 



^/En + M [En -M) = V2n\Ml {N + -) 



(3) 



where N = 2n + i, is the total harmonic oscillator quantum number, n is the radial quantum 
number and I is the orbital angular momentum. Hence the eigenenergies have the same 
degeneracies as the non-relativistic harmonic oscillator. This eigenvalue equation is solved with 
Mathematica, 
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where B{Aiy) 
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— — ^(A^ + §)• The spectrum is non- linear in 
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larmonic oscillator; i.e., the relativistic harmonic oscillator is 



not harmonic. However for small A^ (large M), the binding energy. 



EN-M^M{^ + ...)^hoj{N + ^), 



(5) 



in agreement with the non-relativistic harmonic oscillator. For large A]\f (small M) the spectrum 
goes as 

En^M{AI + ^ + ---), (6) 

which, in lowest order, agrees with the spectrum for Af — )• }13j . 



2.2. U(3) Generators 

The relativistic energy spectrum has the same degeneracies as the non-relativistic spectrum, 
even though the dependence on A^ is different. This suggests that the relativistic harmonic 
oscillator has a higher U(3) symmetry. The non-relativistic U(3) generators are the orbital 

angular momentum £, the quadrupole operator q^ = ■f;jY^J\{2M'^uP'[rr]m + bp]rn ), where 



[rr]m means coupled to angular momentum rank 2 and projection m, and the total oscillator 
quantum number operator, Mnr = J- {2M'^uj'^r'^ + p"^) — |- They form the closed U(3) 
algebra 

[MNRj] = [MNR,qm]=0, (7) 

[li]^'^ = -V2 £6t,i, [/,g]W = -Ve q.m St,2, [q,q]^''^ = 3^/l0 £ 5t,i, (8) 

with Mnr generating a U(l) algebra whose eigenvalues are the total number of quanta N 
and £, qm generating an SU(3) algebra. In the above we use the coupled commutation 
relation between two tenors, t[^^\t!^^^'> of rank ti,t2 which is [T^ ^\r2^*'^]W = [rf^^ T!l^^]{i) _ 

The relativistic orbital angular momentum generators L are given in Eq. ([2]). We shall now 
determine the the quadrupole operator Q^a and monopole operator N that commute with the 
Hamiltonian in Eq.([T]). In order for the quadrupole generator 

n -( {Qm)ll {Qm)l2^-P \ (c>\ 

"^"^ ~ [ a ■ p {Qm)21 a-p{Qrn)22a-p J^ ^^^ 

to commute with the Hamiltonian, [Qm, H] = 0, the matrix elements must satisfy the conditions, 

{Qm)l2 = {Qm)21, (10) 

2[{Qm)luV] + [{Qm)l2,P^]=0, (11) 

2[iQm)l2,V] + [{Qrn)22,P^]=0, (12) 

{Qm)u = {Qm)l2 2(y + M) + {Qm)22 p" ■ (13) 

One solution is 

^--^'\ a-pMu^-Vri^ b# j' ^^ 

where A2 is an overall constant undetermined by the commutation of Qm with the Dirac 
Hamiltonian. 

For this quadrupole operator to form a closed algebra, the commutation with itself must be 
the orbital angular momentum operator as in Eq. ([8]). This commutation relation gives 

y a ■ p £ ^ J 

(15) 



and we get the desired result if A2 = J TT^r^lTTTJ) • The quadrupole operator then becomes 



^ _ / 3 ~ / Ml^2 ^Mu^ ^2 + 2M)[rr](^) + [pp]^^^ Mco^ [rr]^^^ ^'P] hf,\ 

'^"'~^Mu;^h\H + M) [ a-pMio^[rrf} [pp\g^ J^ ' 



which can also be written as 



which reduces to the non-relativistic quadrupole generator for H — )• M. In the original paper 
that derives the quadrupole generators there are two typos. In Eq (6) of that paper [15j, —^ r^ 
should be replaced by Muj"^ r^ and in the non-relativistic quadrupole operator M^a;^[rr]m should 

be replaced by 2M'^uP'[rr]m . Also, the expression for S(yljv) in that paper has a misplaced factor 
of 2 in the denominator. 

For the monopole generator, we can solve the same equations. This has been done [T^. But 
there is a simpler way. From Eq ([3]) we get, 

^^VFTM(£^M)_3^ (18) 

hV2M^ 2 



In the non-relativisitc limit, H + M ^ 2M and the non-relativistic Hamiltonian (H — M) — ?• 
Tujj{N + |) which gives the correct result. 

The commutation relations are then those of the U(3) algebra, 

[M,L] = W,Qm]=Q. (19) 

[L,L]^'^ = -V2 L 5t,i, [L,Qp = -V^ Q 5t,2, [Q,Qp = ^Vw L 5t,i. (20) 

The spin generators in Eq.Q, S, commute with the U(3) generators as well as the Dirac 
Hamiltonian, and so the invariance group is U(3)xSU(2), where the SU(2) is generated by the 
spin generators, [5,5]^*^ = — \/2 S 6t^i. 

3. The Spherically Symmetric Dirac Hamiltonian with Pseudospin Symmetry 

The Dirac Hamiltonian with pseudospin symmetry is [4j 

H = a-p + pM + {l-P)V{r), (21) 

which explains the pseudospin doublets observed in nuclei [Ij. This pseudospin Hamiltonian can 
be obtained from the spin Hamiltonian with a transformation 

75 = ( J l^andM^ -M, (22) 

which gives the pseudospin and pseudo-orbital angular momentum generators jllj 



3.1. The Dirac Hamiltonian with Pseudospin Symmetry and Harmonic Oscillator Potential 
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With the harmonic oscillator potential V{r) = —^ r"^ the eigenvalue equation in the pseudospin 



limit is [To 



Ef;^-M {Ej^ + M) = V2^VM {N + -) (24) 



where A^ = 2n + £, is the pseudo total harmonic oscillator quantum number, h is the pseudo 
radial quantum number and i is the pseudo-orbital angular momentum. While n is the number 



of radial nodes and i the rank of the spherical harmonic of the upper Dirac radial amplitude, n 
is the number of radial nodes and £ the rank of the spherical harmonic of the lower Dirac 
radial amplitude. Again the eigenenergies have the same degeneracy pattern as the non- 
relativistic harmonic oscillator in the spin symmetry limit. This eigenvalue equation is solved 
on Mathematica, 
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where B{An) 
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and A 



N 



V2hLL 

M 



{N + 1). The spectrum is non-linear in 



contrast to the non-relativistic harmonic oscillator; i.e., the relativistic harmonic oscillator is 
not harmonic in either limit. Even for small Aj^ (large M), the binding energy, 
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and hence goes quadratically with the the total pseudo-number of quanta and is non-linear even 
for large M. For large j4^ (small M) the spectrum goes as 



which, in lowest order, agrees with the spectrum for spin symmetry. 



(27) 



3.2. Pseudo-U(3) Generators 

The pseudo-U(3) generators which commute with the Dirac Hamiltonian, [i/jS"] = [H,L] 
[H,Qm\ = [H.,J\f\ = 0, are then obtained by the transformation in Eq ([22|) and are given by 
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The commutation relations are then those of the U(3) algebra. 
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(28) 
(29) 

(30) 
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The pseudospin generators in Eq. ([23|) . S, commute with the U(3) generators as well as 
the Dirac Hamiltonian, and so the invariance group is pseudo-U(3)xpseudo-SU(2), where the 

pseudo-SU(2) is generated by the pseudospin generators, [S,S]^^> = — \/2 S 6t^i. 

4. The Axially Deformed Dirac Hamiltonian with a Harmonic Oscillator Potential 

The relativistic non-spherical harmonic oscillator has also been solved analytically in these two 
limits [To]. The non-relativistic axially symmetric deformed harmonic oscillator will have a 
U(2)xU(l) symmetry. Likewise the relativistic axially symmetric deformed harmonic oscillator 
will have a U(2)xU(l) symmetry in the spin symmetry limit and a pseudo-U(2)x pseudo-U(l) 
symmetry in the pseudospin limit. This will be discussed in a forthcoming paper. 



5. Pseudospin Symmetry Limit and Perturbation Theory 

Even though the harmonic osciUator eigenenergies for the pseudospin symmetry hmit, Eq ()25p . 
are positive, the eigenstates [lOj are not those of Dirac valence states but of Dirac hole states. 
This follows from the fact that the eigenstates have the radial nodal structure of hole states 
rather than valence states. For example, the state with n = i = 0, has by definition a lower 
amplitude with zero radial nodes. This state corresponds to a Ipi state. Using a theorem 

2 

which relates the radial nodes of upper and lower amplitudes [16], the lower amplitude of a Ipi 
state would have to have one radial node if it were a valence state. Hence this state is a hole 
state. This means that perturbation theory can not be used to describe nuclei even though 
nuclei have approximate pseudospin symmetry. In general, the pseudospin limit does not have 
bound valence states independent of the form of V{r) \TU] and hence perturbation theory is not 
possible. For this reason, the analytical solution of the Dirac Hamiltonian with general scalar 
and vector harmonic oscillator potentials would be useful. 

6. Spin Symmetry Limit and Anti-nucleons 

The fact that a nucleon moving in the mean field of a nucleus has a vector and scalar potential 
almost equal in magnitude but opposite in sign has been substantiated by experimental evidence 
in nuclei [1] . An anti- nucleon moving in the mean field of a nucleus will then have a vector and 
scalar potential almost equal in magnitude because charge conjugation leaves the Lorentz scalar 
field invariant while changing the sign of the Lorentz vector field. Thus pseudospin in nuclei 
predicts spin symmetry for an anti- nucleon in a nuclear environment [6l [1] . 

7. Summary 

We have shown that a Dirac Hamiltonian with equal scalar and vector harmonic oscillator 
potentials has an U(3)x SU(2) symmetry and with scalar and vector harmonic oscillator 
potentials equal in magnitude but opposite in sign has a pseudo-U(3) x pseudo-SU(2) symmetry 
and we have derived the corresponding generators for each case. If speculation that an anti- 
nucleon can be bound inside a nucleus is valid [9], the anti-nucleon spectrum will have an 
approximate spin symmetry and, most likely, an approximate U(3) symmetry, because the vector 
and scalar potentials are approximately equal and very large [1]. 
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